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Abstract 

The trends with mass number are examined for the odd-even-staggering (OES) in nuclear binding 
of neutrons and protons through the conventional measures A^ 3 \ The large differences previously 
observed between these trends for even and odd values of these measures is found to arise, in part, 
from the slow variation of binding energies with mass and charge which provides a background 
contribution. This background is estimated with the liquid-drop model, and accounts for the 
greater difference found in proton removal relative to neutron removal. The differences which 
persist after backgrounds are removed can not be treated in the conventional liquid-drop model 
but require the addition of a new term. Such a term is investigated, and its effect on specific values 
of the OES is calculated. The liquid-drop fitting is also applied to a set of separation energies 
constrained to match the specific set of nuclei used to determine the observed values for the odd 
A^ 3 ). The resulting fit for the pairing term is compared to the average value of even and odd 
measures. The effect on this value of the new liquid-drop term is observed, and the change in 
background when the new term is included is also used as an alternate method for determining the 
difference between trends of the even and odd values of the OES. 
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I. INTRODUCTION 



The abrupt changes in binding energy as one goes from a nucleus with an even number 
of neutrons (or protons) to its neighbor with a odd number of such nucleons is know as 
odd-even-stagger (OES). This effect has been much studied and is often associated with 
BCS pairing. However it is known to involve other effects. 1|. 

The experimental observations of OES have been quantitatively examined by several 
measures. One of these is A®, which, for changing neutron number N at fixed Z, is defined 
in terms of binding energies BE by 

A (3) (iV) = ±.5(BE(N + 1) + BE(N - 1) - 2BE(N)) (1) 

An analogous definition /S.^\Z) applies to changing proton number with fixed N. The plus 
sign is used when N is odd, and minus when N is even. This convention guarantees that all 
measures are positive when nuclei with even neutron numbers are more bound than those 
with odd. 

The observed trends with mass number A show a marked difference between cases with 
N odd, A { o\N), and those with N even, Ae (iV). The fits to an A-dependence for neutrons, 
reported in ref.l, were found to be (.82 + 24./v4)MeV and (.94 + 41./A)MeV respectively. 
For protons the fits are found to have an even greater difference, (.75 + 25./v4)Mev and 
(1.11 + 59./A)Mev respectively. If the staggering arose solely from BCS pairing one would 
expect no difference between the behavior for even and odd cases. The purpose of this paper 
is to further examine this difference making use of the liquid-drop model. We treat here the 
situation for both neutron and proton removal. 

As seen in Eq. 1 the measures A^ 3 ) are simply proportional to the second difference with 
respect to neutron number ( or proton number). In addition to the OES, there are also slower 
mass trends involving variations of binding energy with neutron and proton number. These 
variations can also contribute to the second differences. Here, we call the contribution from 
this type of variation "background" . The OES contribution is positive for both even and odd 
cases of A^ 3 ^. The background contribution, however, changes sign due to the convention 
indicated in the definitions. However there would be little difference in its magnitude in 
going from one nucleus to its neighbor if the background variation is slow enough. 

In order to distinguish the specific OES contributions from background effects we label 
these by the short-hand notations ne and no for neutrons and pe and po for protons. For 



the respective background contribution we use the symbol (3 such that 
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When forming the average of the even and odd measures, .5(Ae + A ), the contributions 
from the background f3 cancel. Therefore this average is affected only by the contributions 
from OES. On the other hand, in forming the difference of the measures, e.g., .5(A^ — A®), 
the result is strongly influenced by the background (3. 

II. BACKGROUND 

We use here the liquid-drop model to assist in exploring the differences between the even 
and odd measures. This model is particularly suited for describing the slow variation leading 
to the background. 

The standard liquid-drop model, or more simply the Weizsacker mass formula, was de- 
vised to deal in a simple phenomenological way with physical features of nuclear binding 
energy that are common to classical systems as well. These include a bulk interaction which 
is treated by a volume term proportional to the number of nucleons A. In addition, surface 
tension and Coulomb interactions are treated by two additional terms proportional to A 2 ! 3 
and Z 2 /A 1 / 3 respectively. To account for some of the quantum and isospin effects peculiar 
to nuclei, an asymmetry term is also required. This is taken proportional to (N — Z) 2 /A. 
Another feature of nuclear binding energies is the OES. To a limited extent this feature 
has, until now, been included in the liquid-drop mass formula through a pairing term taken 
proportional to 0.5(( — 1)^ + ( — l) z ) multiplied by function of mass number [i(A). While the 
pairing term of the standard mass model can indeed incorporate the effects of BCS pairing, 
it cannot incorporate some of the features arising from other physics. In particular it cannot 
treat the physics leading to the difference in the mass trends between A® and This 
requires a further enhancement of the liquid-drop mass formula. By studying the observe 
difference in the mass trends of experimental data for odd and even values of A® we can 
also learn about the nature of such an extension of the liquid-drop approach. 
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Because the liquid-drop model provides simple analytic expressions for each of the four 
slowly contributions, it is convenient to evaluate the second partial derivatives of these with 
respect to either neutron number N or charge Z . These can be immediately associated 
with the second differences of A^ 3 ) which contribute to the backgrounds (3 n and (3 p . 

The function p(A) in the pairing term is frequently assumed to have the form c/A 1 ^ 2 . 
This can be replaced by the form p(A) = c\ + C2/A with little affect on the other terms and 
no effect on the rms deviation of the overall fit to the table of measured masses. As pointed 
out in ref.l, the later form is more transparent for indicating the physical origins of this 
term. 

III. NEW EXTENSION OF LIQUID-DROP MODEL 

The liquid-drop model as described above is unable to provide a difference between the 
OES contributions associated with even and odd cases, ne and no, and pe and po, since p 
is assumed to vary little from one nucleus to its neighbor. In order to be able to obtain a 
difference, in the context of the liquid-drop model, an interesting new effect must be added. 

In ref. 1 we examined possible physical sources which lead to such differences. This was 
done from the point of view of mean-field models. To summarizes, we found that a source 
for these differences arose from the degeneracy of the time reversed orbits. This degeneracy 
introduces a non-smooth progression of single particle energies with the successive filling of 
nucleon orbits. That is, the single-particle energies change only after a pair of orbits are 
filled. The interaction terms which reflect the wave-functions of these orbits show a similar 
lack of smoothness with progressive filling. The expected BCS pairing, on the other hand, 
does not provide a difference. 

For convenience we define as 5 n and 5 P the expected non-zero difference between even 
and odd OES, so that 

5 n = ne — no (3) 
5 P = pe — po 

To see how 5 can be incorporated into the liquid-drop treatment we next examine the 
neutron case as a specific example. We add a new smoothly varying part n to the liquid- 
drop formula. This would contributed to the background of ne and no with a magnitude 
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■5d 2 rj n /dA 2 but with opposite signs. 

To link the r) n to ne and no we consider the following expressions for each, involving two 
contributions to the OES. 

ne = no + 5 = (j, - .5d 2 n/dA 2 (4) 
no — fj, + .5d 2 n/dA 2 

These expressions reflect the additional background from n, and /i is the average of ne and 
no. Taking the difference we have 

S = -d 2 n/dA 2 (5) 

This gives a direct link between n and 5, and shows that a non-zero value of 5 requires a 
non-zero function n(A). 

IV. DATA 

We turn next to the measured values of A®. If the full background (3 consists of a part 
Po from the conventional liquid-drop and the additional new term .5d 2 n/dA 2 , we obtain 

A® (iV) = ne - (3 0n - .5d 2 Vn /dA 2 (6) 
Af\N)=no + p 0n + -5d 2 Vn/dA 2 

Af\P) = ne- (3 0p - .5d 2 Vp /dA 2 

Af(iV) = no + (3 0p + .5d 2 r] p /dA 2 

From Eq.2 and the link between r] n and 5 n , Eq.5, we obtain the difference 

.5(Af(A0-Af (A0) = -A,„ + 5n (7) 

(Note that had r] n (A) been missing, the 5 n on the right hand side of Eq. 7 would have been 
multiplied by a factor of 0.5.) Similar expressions follow for protons. 

The experimental mass trends for the left hand side of Eq.7 have been presented in ref.l. 
These trends have been fit by (ci+cz/A) and the explicit values c\ and C2 have been reported, 
having the values of (0.06 + 8./A)MeV for neutrons and (0.18 + 17./A)MeV for protons. 

We estimate here the mass trends for (3 using the following procedure. We fit coefficients 
for the standard 5-term liquid-drop model to the complete table of measured nuclear binding 
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energies For this procedure we used the on-line program [3]. The results for the fitting 
coefficients are 15.7, 17.6. 23.4 and 0.71 for volume, surface, asymmetry, and Coulomb 
terms. These are similar to those commonly reported in text books 4]. 

The volume term, proportional to A, has a vanishing second derivative so to calculate 
Po we need only the remaining three. With the best fit coefficients, the second partial 
derivatives with respect to neutron number of the respective A and Z functions are easily 
found. The values Po n and Po were evaluated for the explicit set of iV and Z used to 
obtain the measured Aa (N) and Aq (P). These sets of values were then fit by a two-term 
expression c\ + Cij A giving the resulting trends for /3 

p 0n = (-0.0066 - 16A3/A)MeV (8) 
P 0p = (-0.087 - 32.79/ 'A)MeV 

From Eq.7 we then estimate the mass trends for 8 n and 5 P using the coefficients provided 
in table 3 of ref.l for the observed difference of measures, 

S n = (0.05 - 8./A)MeV (9) 
5 P = (0.09 - l6./A)MeV 

The full background Pq — .55 then gives 

p n = (-0.04 - 12./ A)MeV (10) 
P P = (-0.14 - 25. /A)MeV 

The greater values for the proton background result in the greater difference between the 
observed even and odd measures. 

Using this dependence for P and the values of 5 we deduce 

(11) 

no = (0.82 + 24. /A) + (0.04 + 12. /A) = (0.86 + 36. /A) 
ne = (0.86 + 36. /A) + (0.05 - 8. /A) = (0.91 + 28. /A) 
po = (0.75 + 25. /A) + (0.14 + 25. /A) = (0.89 + 50. /A) 
pe = (0.89 + 50. /A) + (0.09 - 16. /A) = (0.98 + 34./A) 

Finally, as a numerical check, the average of ne and no and the average of pe and po give the 
same results as the average of the respective Ae and Ao as they must, by construction. 



V. SEPARATION ENERGIES FOR CONSTRAINED SETS OF NUCLEI 



The measure A( 3 ) is also proportional to the difference between neighboring separation 
energies. Therefore we studied fits to the sets of measured neutron and proton separation 
energies. For this we also used the liquid-drop approach but fit the coefficients of the 
various terms to the sets of respective separation energies (neutron and proton) rather than 
the binding energies. In this study we have included only nuclei which are actually involved 
in obtaining the sets of values for A (N) and A (P) reported in ref.l. The constraints on 
the data are outlined there and, corresponding to these, for neutrons we have used separation 
energies for nuclei with Z even, N > Z, and Z greater or equal to 8 and N greater or equal 
to 11. For protons we have used N even, N > Z, N great or equal to 12 and Z greater or 
equal to 9. Because of the close link between OES and the separation energies this promised 
to provide additional confirmation of the approach to backgrounds taken above with the fit 
to binding energies. 

With the constrained set of nuclei the rms deviations for the best fit to the separation 
energies, with the 5-term conventional form of the liquid-drop formula, was 0.62 MeV for 
neutrons and .57 MeV for protons. (Note, this is the rms from the best fit to the separation 
energies and not binding energies.) The coefficients of the pairing term was 12.4/A 1//2 MeV 
for neutrons and 13.9/A 1 / 2 MeV for protons. When the form of the pairing term was taken as 
(ci+cz/A) the best fit was given by (0.897+29.7/A)MeV for neutrons and (.96+37.2/A)MeV 
for protons. These fits provide the same rms deviations. In Section III we found the average 
value of the odd and even measures for neutrons and protons to be was (0.88 + 33./A)MeV 
and (0.93 + 42. /A) MeV as reported found in Ref.l. This is in relatively close agreement 
with this direct fit to the sets of separation energies using only the four slowly varying terms 
of the liquid-drop. 

Since the 5-term fit lacked the flexibility to provide a difference between no and ne, for 
example, it could be expected that the rms deviations might be improved by repairing this 
deficiency with the additional term rj discussed above. One of the contributions to 5 is 
anticipated to be proportional to the level spacing of the single-particle energies. A good 
fit to the mass trends of this quantity is associated with the form 1/A. Consequently one 
might anticipate 77(^4) to have the form Aln(A) in order that its second derivative might 
provide S. 
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When a additional term of the form (cAln(A)) was added to the 5-term liquid-drop 
formula, the rms deviation value was indeed reduced to .58 MeV for neutrons and .54 MeV 
for protons. The values of the best fit coefficients for the additional terms are found to be 
c = —9.98 and c = +18.75 for neutrons, and protons respectively. However the coefficients 
of the 4 traditional terms change from their original values when the extra r] terms were 
added. 

We next considered the values of S n and 5 P which might be obtained from d 2 rj/dA 2 . If 
the added terms (cAln(A)) did not affect the other coefficients then the fitting coefficients 
themselves would lead directly to the form of 5. However this is not the case, and the values 
of /3q are perturbed in the new fitting process. To incorporate this change we consider the 
"effective" value of .5d 2 /dA 2 to be the difference between the background for the augmented 
liquid-drop (with the rj added) the traditional liquid-drop (without 77). We thus consider 

.5d 2 r]/dA 2 = A(3 Q + 0.5c/ A (12) 

The net change A(3 when 77 is added is (-0.038 + 9.72/A)MeV and (-.04 - 4.73/A)MeV 
for neutrons and protons respectively, and the constants c are coefficients of Aln(A) from 
the fits. The value of 5 is then given by multiply .5d 2 i]/dA 2 by -2. 

This procedures produces values for 5 of (0.08 - 9.5/A)MeV and (0.08 - 9.3/A)MeV 
respectively for neurons and protons. This is qualitatively similar to the differences 5 (0.05 — 
8./A)MeV and (0.09 - 16./A)MeV which were estimated in Eq.ll of Section IV. The l/A 
terms in the expressions for 5 appear to have negative coefficients. If the physics behind this 
term is, as suggested in Ref. 1, the sum of a single-particle level spacing term (providing a 
positive coefficient of order 45) and potential energy term (providing a negative coefficient), 
then the potential energy term seems to dominate but be of similar magnitude as the level 
spacing term. 

The addition of the Aln(A) also leads to direct fits for the pairing function which are 
slightly closer to the observed average of A^ and A^ reported reported in Ref.l. For 
neutrons the fit provided (.88 + 33./A)MeV, precisely the observed value. For protons it 
provided (.96 + 38./A)MeV, slightly closer to observed value than without the additional 
term. 
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VI. LIQUID-DROP PAIRING TERM FITS 



We also examined the sensitivity to the pairing term when the liquid-drop formula is 
applied to the set of separation energies values used in forming A^\N). For this purpose 
we performed a series of least-square fits starting with only the pairing term ( and none of 
the slowly varying terms) and then added, one-by-one, the other four terms of the model. 
When all 5-terms of the liquid-drop are used the fit to the pairing term provides a coefficient 
of about 12.4/A 1//2 quoted above. When no terms except the pairing term is used the fit 
gives a value about 13. /A 1 / 2 with a huge rms deviation. As more and more of the terms are 
added the fits to the pairing term trend slowly downward toward the fit obtained with the 
full complement of terms, and the rms deviations decrease substantially. What is remarkable 
is the small range over which the best fit pairing coefficient changes. This suggests that the 
pairing term is strongly determined by the set of separation energies with little influence 
from the other important terms. Those terms are, of course, quite important for obtaining 
a reasonable fit. In Table I we demonstrate this by providing the pairing coefficients and 
the rms deviations for fits involving increasing numbers of liquid-drop terms. 

TABLE I: Fits of c for the coefficient of the pairing term c/A 1 / 2 with increasing number of terms 
in liquid-drop formula. The fit is to neutron separation energies for even Z, N > Z, Z greater or 
equal to 8, and for two values of N m i n , 10 and 11 



Terms used 




10 




11 


rms 












MeV 


pair. 


12.9 




13.1 




8.02 


pair.+vol. 


12.7 




12.7 




1.81 


pair.+vol.+surf. 


12.7 




12.7 




1.81 


pair . +vol.+surf. +coul. 


12.6 




12.5 




1.39 


pair ,+ vol . +surf. +coul .+ asym 


12.4 




12.4 




0.62 



We next examined the correspondence between fits in which the pairing coefficients has a 
2-term form (ci + C2/A) and those in which it has the 1-term form (c/A 1 / 2 ). We considered 
the set of 2-term pairing coefficients (ci + c 2 /A) for which the values of A were taken for 
each nucleus belonging to the set of nuclei used for forming the A^\N). We then fit this 
set of values with a curve employing the 1-term form c/A 1 / 2 . The results of this exercise for 
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values of c giving the best fit are show in Table II. 

TABLE II: Fit by the curve c/A l l 2 to values of c\ + 02/ A for each of nuclei associated with (N). 
The values shown are c for each combination of c\ and c 2 



value of ci 


c 2 =29. 


c 2 =30 


c 2 =31 


0.89 
0.90 


12.24 
12.34 


12.35 
12.45 


12.46 
12.56 



This seems to explain why similar fits result with the use of either the 2-term or 1-term 
forms. It also suggests that a change of 0.1 in the value of the coefficient of the 1-term fit 
corresponds roughly to a change of 1. in the value or c 2 or .01 in the value of c\. 

VII. SUMMARY AND CONCLUSION 

The mass trends of OES have been examined using the measures A®. The observed 
binding energies indicate a large difference been even and odd values of these quantities. It 
is pointed out that the slow variations of the binding energies, with mass and charge, can 
significantly influence this difference by contributing background contributions f3 in addition 
to the OES. We have studied this influence by using the results of the liquid-drop model to 
describe these variations. While the backgrounds are significant for all the measures, they 
are found to be greater for the proton removal than for neutron removal. Consequently this 
can account to the greater differences in the observed values. This effect probably arises from 
the greater influence of the Coulomb energy and the asymmetry energy with a predominance 
of neutron number over proton number in the mass tables. 

After the background is removed differences in OES are still anticipated from effects 
discussed in ref. 1, and such differences 8 cannot be treated by the conventional terms 
employed in the liquid-drop model. Consequently we propose adding to the liquid-model a 
new contribution arising from the specific physics contributing to this difference between even 
and odd trends in the staggering. We have examined the new contributing function, labeled 
here as r)(A), and find that d 2 r]/dA 2 = —5. We have combined the observed difference 
.5(Ae — Ao ) and the liquid-drop background to estimate 8 for neutrons and protons. 
These estimates provide (0.05 — 8./A)MeV for neutrons and (.09 — 16./A)Mev for protons. 
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For comparison, the difference between Ae and A provides the much greater values of 
(0.12 + 16./y4)MeV for neutrons and (0.36 + 34./A)MeV for protons. Note the sign change 
in the coefficient to the 1/A terms. 

We have also applied the liquid-drop model fit directly to the separation energies for 
the specific set of nuclei used to constructed the observed A ■ This procedure, without 
the additional contribution 77 led to a pairing function close to the observed averages of 
the odd and even OES mass trends. We have also added a contribution 77, assuming a 
form Aln(A) motivated by the expectation that 5, and hence d 2 rj/dA 2 , might follow a form 
1/A. This addition improved the rms fit and also brought the pairing term function closer 
to the observed average of even and odd OES; it providing an exact match for neutron 
removal and qualitative match for protons. Adding the 77 term also led to a change of the 
background contribution which could be associated with an effective value of d 2 r]/dA 2 , and 
consequently 5. This method of determining 5 provided values in qualitative agreement with 
those obtained by removing an estimated background: (0.08 — 9. /A) compared to (0.05-8. /A) 
for neutrons, and (0.08 — 9. /A) compared to (0.09 — 16. /A) for protons. 

A study of the liquid-drop fit to the observed separation energies suggests that the data 
strongly determines the liquid-drop pairing term values, and that these values are little 
influenced by the presence or absence of the other terms in the model. Namely, the best fit 
for the pairing term changes from 13.0/ A to 12.4/ A as the number or other liquid-drop terms 
ranges from zero to the conventional four. We also made a detailed comparison of the two 
forms for the pairing function, c/A and (ci + C2/A), where a change of 0.1 in c corresponds 
to a change of 0.01 in c\ or 1.0 in c^. 
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